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CONVERTING FROM 3.6 AND 4.5 FLUXES TO STELLAR MASS 
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ABSTRACT 

We use high spatial resolution maps of stellar mass and infrared flux of the Large Magellanic 
Cloud (LMC) to calibrate a conversion between 3.6 and 4.5 /im fluxes and stellar mass, Af* = 



10 5.65 F 2.85 



FZ-PFrh 86 (tMs) M , where fluxes are in Jy and D is the luminosity distance to the source 



4-5 V 0.0.5/ 

in Mpc, and to provide an approximate empirical estimate of the fractional internal uncertainty in 
M* of 0.3y/N/10 6 , where N is the number of stars in the region. We find evidence that young stars 
and hot dust contaminate the measurements, but attempts to remove this contamination using data 
that is far superior than what is generally available for unresolved galaxies resulted in marginal gains 
in accuracy. The scatter among mass estimates for regions in the LMC is comparable to that found 
by previous investigators when modeling composite populations, and so we conclude that our simple 
conversion is as precise as possible for the data and models currently available. Our results allow for 
a reasonably bottom-heavy initial mass function, such as Salpeter or heavier, and moderately disfavor 
lighter versions such as a diet-Salpeter or Chabrier initial mass function. 
Subject headings: galaxies: fundamental parameters, Magellanic Clouds, stellar content 



1. INTRODUCTION 

The measurement of the mass of the stellar population 
in a galaxy, M*, is a direct, although integrated, mea- 
sure of the star formation history of the system. We are 
increasingly realizing that is intricately connected to 
other galaxy properti es, such as current star formation 
rate and morphology (Kauffman n et al.l l2~004). Unfortu- 
nately, the measurement of M„ is indirect and subject to 
significant systematic uncertainties. 

There are effectively two principal approaches to mea- 
sure M*. First, one can measure the dynamica l mass of 
a ga laxy, via kinematic s ([Cappellari et aLll2006D or lens- 
ing (jAuger et al.ll2006T) . and then somehow model and 
subtract the contribution of dark matter to that mea- 
sured mass. This approach is predicated on the success- 
ful subtraction of a l arge unseen mass, which i s dominant 
in galaxies overall ([Zaritsky fc White] 1994) an d even 
substantial within the optical radius ([Cappellari et al.l 
2006| I Auger et al J 120061 ). and can easily lead to uncer- 
tainties of the order of the measurement itself. Sec- 
ond, one can rel y on stellar population models (for 
example those of iBruzual &; Charlotl 120031 ) to connect 
M* to an observable, such as the luminosity in a se- 
lected passband, a color, or the spectral energy distri- 
bution as obtained either from spectroscopy or multi- 
band photometric observations. In this case, one is 
also subject to modeling uncertainties, although here 
those arise from uncertainties in the stellar initial mass 
function (IMF), the stellar evolution models, partic- 
ularly in the stages where stars are at their most 
luminous (for some examples drawn from an exten- 
sive literature on the topic see [Langer fc Maederl ll99E 



Maraston et al. 2006; Conroy et al. 2009; McQuinn et al. 
2011tlDalcanton et al.ll20l"li r and the star formation his- 
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tory. Careful use of scaling relations (|Bell et al.l 120041) . 
which are themselves indirect measurements of mass, 
provide a bridge across the two methods, greatly enlarge 
the sample size and help mitigate these problems, but 
result in prescriptions that may be accurate on average 
but highly uncertain in specific cases. 

Having these two independent approaches to measure 
M* is advantageous in that the results can be com- 
pared to uncover any systematic diffe rences. For example 
McLaughlin & van der Marcl (2005) compare the results 
obtained for stellar clusters and conclude that there are 
no egregious differences. The magnitude of the differ- 
ences are now constrained to be relatively small (few 
tens of percent) and therefore modeling subtleties are 
critical and difficult to control for. We address this is- 
sue by presenting and applying a third method, which 
is more closely related to the photometric method out- 
lined above but avoids several perilous assumptions. Our 
aim is to calibrate an easily observed quantity that is 
available for large numbers of galaxies and relatively im- 
pervious to the effects of extinction, ongoing star forma- 
tion, and details of rare, but luminous phases of stellar 
evolution, and to determine the uncertainties introduced 
by those phenomenon. We choose here the more mod- 
est goal of identifying and calibrating such an estimator 
of M„ due to our interest in applying such an estima- 
tor, rather than the more demanding goal of fully un- 
derstanding the physical origin of the scatter and any 
possible systematic biases we uncover. 

To help us measure M*, we select the 3.6 and 4.5 
/im fluxes, F3.6 and F4.5, which are beco ming increas- 
ingly available for large samples of g alaxies dSheth et all 
120101 ) with t he advent of the Spi tzer ([Werner et alJl2004f > 
and WISE ([Wright et all [2010) telescopes, and which 
are nearly minimally sensitive to young stellar popula- 
tions and dust absorption and emission. We will test 
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combinations of -F3.6 and -F4.5 as proxies for M* in one 
galaxy, the Large Magellanic Cloud (LMC), whose ex- 
tinction and sta r formation hi s tory have been mapped 
in spatial detail (jZaritskvl 11999: Harr is fc Zaritskvll2009f ) 
and that has been extensively observed with the Spitzer 
telescope (jMeixner et al.l 120061 ) . 

The basic concept we exploit is that we can use 
the star formation histories (SFHs) recovered from syn- 
thesizing the stellar optical color-magnitude diagrams 
(CMDs) that are available for nearby galaxies ( e.g. 
lHarris fc Zaritskvl l200l 120091 : iDalcanton et aD 120091 ) to 
calculate the stellar mass on a region-by-region basis and 
use those measurements to calibrate F3.6 and .F4.5 as 
tools with which to measure M* . We then use the scatter 
in that correspondence to uncover any additional param- 
eters that can be used to refine the measurement and to 
determine the underlying uncertainty in this measure- 
ment. Although this approach shares some of the diffi- 
culties faced by the stellar population synthesis approach 
(dependence on initial mass function and on stellar evo- 
lution models), it eliminates the poorest constrained as- 
pect of those models, the star formation history, and 
minimizes the effect of poorly understood rare phases 
of stellar evolution, such as the thermally pulsing AGB 
phase experienced by intermediate mass stars and the 
red, core H e burning phase experien ced by stars of masses 
> 3.5M (|Melbourne et al.l l2012h . The latter advan- 
tage is realized because the CMD modeling depends on 
the number of such stars rather than their luminosities. 
These stars are relatively rare, and hence have little ef- 
fect on the CMD modeling, but extremely luminous, and 
hence have a large effect on global colors and luminosi- 
ties. Recently available SFH maps, as described in §2.1, 
provide the necessary data to calculate resolved stellar 
mass maps, which are then compared to local measures 
01 F3.6 (§2-2), to calibrate i^.g and F4.5 as stellar mass 
tracers (§3). We summarize our findings in §4. 

2. THE INPUT DATA 

2.1. Stellar Masses 

The spatially-resolved SFH of the Large Magellanic 
Cloud ((Harris fc Zaritskvl 12001 hereafter HZ), provides 
the information necessary to construct a resolved stellar 
mass map of the galaxy for an adopted IMF. Because 
the regions are rectangular and the available algorithms 
m to measure a luminosity work with circular 

apertures, we simply use the inscribed circle within each 
of the HZ regions to define the zones in which we will 
also measure the infrared luminosities, and correct for 
the differences in areas. The stellar mass we calculate 
is the integral over mass for the adopted IMF, Salpeter 
(|Salpeterl ll955). normalized by the star formation rate 
at a particular time, integrated over time for the life- 
time of the LMC. We will discuss later (§3) a correction 
for the mass returned by evolved stars to the interstellar 
medium. 

One technical point in this procedure regards the du- 
ration of the earliest (ol dest) bin. Our analysis method 
(Harris & Zaritskv 2001) finds the number of stars in the 

1 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which arc operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



color magnitude diagram that match a population drawn 
from a particular isochrone, and then calculates a star 
formation rate by dividing that number by the length 
of time represented by that isochrone. That duration 
is easier to define when a particular isochrone is brack- 
eted by other isochrones, but is not well defined for the 
oldest isochrone. As such, the star formation rate at ear- 
liest times will depend on the assumption of when the 
LMC started forming stars in this first age bin (e.g. did 
the LMC start forming stars 11 Gyr or 13 Gyr ago?). 
While this issue leads to some uncertainty in the de- 
rived star formation rate at these times (for example, see 
lEskew fc Zaritskvl 1201 If ) it does not pose a problem for 
the current analysis because we integrate over the length 
of the bin (we are interested here in the total number of 
stars rather than in the rate of star formation) . 

The uncertainties in the recovered star formation 
rates from CMD synthe s is are d if ficult to ascertain 
(lHarris fc Zaritskvl 120011 iDolphinl 12001 iWeisz et al.l 
120111 ). Much of that uncertainty is related to how pop- 
ulations are partitioned among adjacent bins in age or 
metallicity. These errors are unlikely to affect our re- 
sults significantly because they arise from the fact that 
such populations are difficult to disentangle using broad- 
band photometry, and hence are unlikely to have highly 
differential properties in our photometry. Uncertainties 
over long times appear to be well controlled, as exem- 
plified by tests done by Harris & Zaritskv (2001) using 
globular clusters, the global patterns in the star forma- 
tion history seen in the LMC, and the correspondence 
of certain feature s in the star formation historie s of the 
LMC and SMC (lHarris fc Zaritskvl l200l [2009h . Such 
errors, even if they are significant, are likely to add ran- 
domly to the relationship between M* and -F3.6 because 
there is little connection between stellar density and star 
formation history in the Clouds, outside of poss ibly in 
the LMC bar region (see lHarris 

jLZaritsky 20Ql). 
The key systematic uncertainty in our approach comes 
from the selection of the IMF. The behavior of the 
IMF at the high mass end [M > few M ) impacts 
th e derivation of the SFH , as discussed for example 
bv lHarris" fc Zaritskv (2001). In determining total stel- 
lar mass, the low mass form of the IMF is a poten- 
tially larger source of uncertainty. Here we assume 
a Salpeter IMF and note that scaling corrections cor- 
responding to other choices of IMF can be calculated 
using population synthesis cod es, such as PEGASE 
(|Fioc fc Rocca-Volmerangelfl997f) . 

2.2. IR Fluxes 

We use the published, calibrated mosa iced images of 
the L MC produced by the SAGE s urvey (jMeixner et alJ 
120061 ) using the IRAC instrument (|Fazio et al.ll2004f ) on 
the Spitzer telescope to measure the fluxes, at both 3.6 
and 4.5 /im, in circular apertures that are matched to 
the regions for which HZ present star formation histo- 
ries. The mosaics have pixels that are 2"on a side, or 
alternatively 9.4016 x 10~ n steradians. The units of the 
surface brightness measurements are MJy/sr. We use 
the IRAF task PHOT, with a background value set to 
zero to calculate the flux in the region. The majority of 
the regions are quite luminous, so the uncertainties are 
dominated by Poisson statistics and much smaller than 
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Fig. 1. — The distribution of Fg.a measurements. The sharp drop 
in distribution near a value of 0.04, shown by the right, dotted, 
vertical line, indicates our upper limit on a possible uniform sky 
value (see text for discussion). Our preferred sky value of 0.02, 
shown by the left vertical line, falls close to values obtained for the 
darkest regions available. 

those associated with the calculation of the stellar mass. 
After integrating over an aperture and converting units, 
we then present fluxes in Jy. 

As mentioned previously, we have artificially set the 
background to zero. Now we attempt to estimate the 
background directly from our results. Our challenge is 
that the LMC extends beyond the region covered by the 
mosaic, so there are no "empty" regions within the im- 
age. First, we examine the distribution of measured 
mean surface brightnesses within the various apertures 
(Figure [1]). We do this rather than examine the image 
directly, searching for dark patches, because we are not 
attempting to estimate the blank field background but 
rather the mean contribution from the blank field plus 
galactic stars plus any galactic diffuse emission. As such, 
selecting the darkest regions could result in a background 
estimate that is biased low. From Figure [TJ we conclude 
that a robust upper limit to the mean background is 
likely to be 0.04 MJy/sr, at which point the number of 
apertures drops precipitously, signaling the lower limit 
of LMC-related emission. A slight peak in values at 0.02 
MJy/sr suggests an alternate choice for the correct back- 
ground level. 

In Figure [51 where we plot the relationship between 
-F3.6 and M*, we show the results of a different approach 
to estimating the background. We calculate the vertical 
scatter about the best fit linear relationship for different 
adopted background values and identify the background 
level that minimizes the scatter. The horizontal scatter 
introduced by measurement errors is negligible given the 
large fluxes over these regions. We find that the best 
fit value is 0.02 MJy/sr. In the three panels, we show 
the relationship for adopted background values of 0.00, 
0.02, and 0.04 MJy/sr. If no background correction is 
applied, then there is a slight deviation of the lower end 
of the distribution away from the 1:1 line. The lack of 
a strong turnover in the relationship demonstrates that 
the background term is not highly significant. When we 
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Fig. 2. — An alternative justification for our choice of background 
level. The left panel shows the relationship between F3.6 and M t 
when no background correction is applied. The middle panel shows 
the result with our preferred background value of 0.02 MJy/sr. 
The right panel shows the result when adopting the upper limit 
of 0.04 MJy/sr. When no correction is applied, the lower part of 
the distribution bends slightly downward below the best fit line. 
When our upper limit correction is applied, we see a tail of points 
that suggests an unphysical, asymmetric increase of the scatter for 
regions with lower total stellar mass. Our preferred value minimizes 
the veritcal scatter about the best fit line (the horizontal scatter 
due to measurement errors is negligible given the large count values 
in these regions). 

adopt our upper limit of 0.04 MJy/sr, we oversubtract 
the background contribution in many fields leading to the 
leftward tail of regions with significant stellar mass but 
little IR emission. For our preferred value, 0.02, neither 
of these problems occurs. A few regions asymmetrically 
scatter to the right in the Figure and there is a systematic 
deviation from the relationship at large flux values, but 
we will discuss the origin of those issues below. We do 
the same analysis for F4.5 and find a best-fit background 
of 0.017 MJy/sr. These backgrounds are subtracted from 
all fluxes discussed subsequently. 

3. RESULTS AND DISCUSSION 

We present the results of our analysis of the fluxes and 
stellar masses for each LMC region in Figured] The up- 
per left panel in the Figure shows the raw result, using 
our adopted background value, for the correlation be- 
tween ^3.6 and M*. The solid line represents the mean 
linear relationship (slope = 1 in our log- log plots). The 
dispersion in values about this line is 0.121 dex, when we 
exclude points that deviate from the line by > 0.3 dex 
(shown by the dotted red lines), and the distribution of 
residuals is reasonably well approximated by a Gaussian. 
This degree of scatter corresponds to mass estimates that 
have la uncertainties of ~ 30%. 

This scatter can, in principle, be reduced if it arises 
from physical sources, and if we can identify and account 
for those sources. One such potential source is the varia- 
tion in star formation history from region to region across 
the LMC. As such, one might suspect that using an inde- 
pendent measure of SFH variations, or mean stellar age, 
to refine the relationship between i^.g and M* would lead 
to lower scatter. We investigate utilizing the one simple 
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Fig. 3. — The dependence of the relationship between -F3.6 and 
M t on IR color. The upper left panel reproduces the relation 
shown in the middle panel of Figure [2] The upper right panel 
shows the residuals from the best fit line in the upper left panel, 
8, as a function of IR color (slope —1.85) The lower left panel 
shows the relationship between F3.6 and M* once we correct for 
the relationship in the residuals with color. The dotted lines are 
±0.3 from the best fit and represent the region used to evaluate 
the dispersions quoted in the text. 

directly observable measurement we have of such popu- 
lation variations, color, in Figure [3J In the upper right 
panel, we show the residuals in the original relationship, 
S, vs. color, \og(F4, 5 /F 3 , 6 ). A linear relationship with 
slope —1.85 is shown in red and suggests that there is 
indeed a weak connection between the residuals and this 
measure of the stellar populations. When we account for 
this correlation between IR color and residual, there is 
modest improvement (lower left panel Figure and a 
dispersion about the best fit line of 0.117, similar to be- 
fore the correction, but the systematic deviation of the 
upper end of the distribution relative to the lower end is 
lessened. 

Given the data we have for the LMC, which we will 
not have in general, we can probe somewhat further into 
the role of stellar populations in affecting F3.6 as a mass 
tracer. Although we will not be able to apply this knowl- 
edge in general, the exercise might help highlight the 
underlying cause for the scatter and provide guidance 
for other approaches aimed at dealing with variations in 
M/L3.6 (|Meidt et all 12012ft . In Figure H we show the 
correlations between the residual, S, and the fraction of 
the stellar population, /, that is younger than 1.3 Gyr 
(upper left panel) and 0.33 Gyr (upper right panel). In 
the lower panels of the Figure, we show the relationship 
between those same fractions and the IR color, which we 
used previously to correct the relationship. The correla- 
tions between 5 and either /1.3 or /0.33 are both strong, 
with the stronger being with the youngest stellar frac- 
tion (Spearman rank correlation coefficients of —0.328 
and —0.493 respectively). Both of these correlation coef- 
ficients correspond to highly significant correlations, with 
probabilities of occurring randomly in a sample of this 
size < 6 x 10 -31 , and they are both stronger than the cor- 
relations seen with color (correlation coefficients of 0.177 
and 0.118, respectively). Interestingly, the residual from 



Fig. 4. — Correlations between the residuals from the mean re- 
lation between F3.6 and M*, <5, and the fraction of stars younger 
than 1.3 Gyr (upper left) or 0.33 Gyr (upper right), and between 
IR color and these same fractions (lower panels). All panels contain 
statistically significant correlations, but the strongest is between 8 
and /0.33 (see text for correlation values). 

the mean trend of -F3.6 vs. Af* is a better diagnostic 
of the fraction of very young stars than is the IR color. 
The finding that /0.33 is more strongly correlated with S 
than color is with S (—0.493 vs. —0.348) demonstrates 
that IR color alone cannot act to correct the relationship 
entirely for stellar population variations. Importantly, 
even relatively modest fractions of young stars (~5%) 
are sufficient to significantly affect M^/F^.q. 

Although we have found that the residuals in the 
i^3.6 — M* relationship correlate with /0.33 and that IR 
color cannot fully correct for this trend, we find that 
correcting for the trend does not significantly lower the 
scatter in the F3.6 — M* relationship. We correct for the 
correlation between S and /0.33 by fitting a linear rela- 
tionship between 5 and log(/o.33) to produce the results 
shown in Figure [5l The corrected data have a scatter 
of 0.115 (30% in mass). Some improvement is visible, 
particularly in the S distributions (lower panels), but the 
decrease in the scatter is marginal. We conclude that us- 
ing even superior data than a single IR color to remove 
stellar population variations (at the level of a linear fit 
to the residuals) is insufficient to substantially lower the 
scatter. Furthermore, this level of scatter is compara- 
ble to that inherent to more involved and data-intensive 
treatments , such as the two-co lor method at H-band de- 
veloped bv lZibetti et al.1 (|2009f ). 

In addition to the issue of scatter about the mean rela- 
tionship, there is a population of regions that consistently 
lie significantly rightward of the mean relation between 
M* and -F3.6 (Figures [3] and [5]). To understand the origin 
of this discrepancy, we plot the location of regions within 
the LMC with large residuals (5 < —0.3) on the 8 /.tm 
imag e of the LMC (Figur e Et, also from the SAGE sur- 
vey (jMeixner et alJ 120061 ). We find that many of these 
outliers correlate with regions of hot dust, as traced by 
the 8 /jm image. However, this correspondence is far 
from ideal and some outliers lie in unremarkable regions, 
highlighting the difficulty in using any color, even one 



5 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 
log F„ [Jy] log Lja [Jy] 




-0.5 0.5 -0.5 0.5 

6 6 



Fig. 5. — The dependence of the relationship between -F3.6 and 
M t on the fraction of stars younger than 0.33 Gyr, /o.33- In the 
left panels we plot the uncorrected relationship between .F3.6 and 
M„ and the distribution of residuals, 8 about that relationship. In 
the right panels we plot the same, except here we have applied a 
best fit correction between 6 and /0.33. Although improvement in 
the relation is evident, it is relatively modest even though we have 
a direct measurement of the fraction of young stars. 




Fig. 6. — Spatial distribution of regions with large residual from 
the -F3.6 vs. Mi, relationship, superposed on the 8 11m mosai c 
image of the LMC from the SAGE survey (Mcixner et aL 2006). 
Specifically, we selected regions with residuals < —0.3 from the 
age-corrected relationship shown in Figure [3] Although many of 
the outliers correspond to regions of high 8 /im flux, about a third 
are in regions that appear to be unremarkable. 

that identifies regions with warm dust, to correct F 3 . 6 
for contamination. The contamination probably comes 
from various components, ranging from wide-scale hot 
dust that is easily visible in the image to compact sources 
such as extreme AGB stars, which can have J — 3.6 col- 
ors m excess of 5 (Bl um et al.ll2006f ) . While the latter are 
evident when one has photo metry of individua l stars, as 
can be obtained in the LMC IjBlum et al.ll2006f) . they will 
be difficult to pick up in integrated populations, in data 
with limited color baselines, and in data with limited S/N 



over small spatial scales. Even with sophisticated t reat- 
ments, such as that described bv lMeidt et al.l (|2012t ). one 
might at best be able to exclude all of these outliers from 
the analysis, at which point one would return to our es- 
timate of the scatter of ~ 0.12 (or 30% in mass for the 
regions considered here). 

One aspect that we have neglected so far is the re- 
turn of mass to the interstellar medium by evolved stars. 
What we have calculated is the sum of the mass of all 
stars formed, which is different than the sum of the 
mass of all current stars. To estimate the systematic 
difference between the two, we use PEGASE models 
(|Fioc fc Rocca-Volmerangelll997f) with the global SFH 
of the LMC as we did in lEskew fe Zaritskvl (poll. We 
find that 94% of all the mass formed into stars by the 
current time remains in stars and remnants. We will 
therefore correct our calculated stellar masses by 0.94. 
In principle, this correction should vary with the SFH of 
each region, but the differences are well below the 30% 
region-to-region scatter, and so we neglect them. 

The color-corrected relationship shown in Figure[3]sug- 
gests that the remaining intrinsic uncertainty, excluding 
the asymmetric tail, is about 0.12 dex, or about 30% for 
regions containing somewhere in the vicinity of 10 5 5 to 
10 7 M Q . If this scatter is random, then estimating the 
masses for larger regions should result in lower fractional 
uncertainty. We can estimate how well these fluctuations 
average out by combining all of our data into an estimate 
of M* for the entire LMC and comparing to the extrap- 
olation of our 1:1 calibrated relationship. To do this, we 
exclude regions that have 5 < —0.3. We find that result 
is exceedingly close to the extrapolated relation (using 
only F3.6 we infer M* = 1.72 x 10 9 M Q and we calculate 
M* = 1.81 x 10 9 M Q using the stellar masses for 1131 out 
of 1180 regions and the correction for mass returned to 
the ISM). If we do not remove any outliers, we obtain 
1.90 x 10 9 and 1.87 x 10 9 M Q , respectively, and similar 
results (1.81 x 10 9 and 1.93 x 10 9 M Q , respectively) if 
we use the relationship with both -F3.6 and -F4.5. Either 
way, this analysis implies that the internal uncertainty 
drops to a few percent when considering systems with 
M»~2x 10 9 M Q . This gain of a factor of at least 10 
in precision corresponds well to the naive \/N gain ex- 
pected in going from regions with ~ 10 6 to ~ 10 9 M Q . We 

conclude that scaling the uncertainty by yN provides a 
reasonable estimate of the expected gain (or loss) in pre- 
cision and suggest that on galaxy scales we are limited 
by systematic rather than statistical uncertainties. 

One final check of our results is possible by comparing 
the dynamical mass, Mjj, of the LMC out to the radius 
sampled here and the stellar mass calculated ab ove. For 
a rotation speed of 87 km s _1 (jOlsen "eTaHIMll) out to a 
radius of 3 kpc, corresponding to the area covered by our 
survey, we estimate that the enclosed mass is 5.0 x 10 9 
Mq . Given the unknown contributions of gas and dark 
matter, plus the o(l) geometric corrections necessary in 
going from our assumed spherical geometry in using M = 
rv 2 to a more realistic disk potential, the two estimates 
are not in significant conflict, and certainly satisfy the 
basic constraint that M* < Mn. 

Our preferred calibrations at the distance of the LMC 
are therefore 

M* = 10 5 ' 97 F 3 . 6 
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when _p4.5 is unavailable and 



AL = 10 



5.65 7^2.85 



^3.6 '4.5 



1.85 



otherwise, where F 3 6 and are fluxes in Jy and M* 
is given in solar masses. The expression for a source at 
arbitrary distance, £>, in Mpc is 



D 

(105 



if we adopt a distance of 50 kpc for the LMC. 

Previous work combining stellar mass esti mates from 
integ rated spectroscopy and photometry (Zh u et ahl 
l2010h provided an alternative calibration of the relation- 
ship between F3.6 and M*. They present two transfor- 
mations, one that depends solely on L3.6 and another 
that also utilizes g — r. We can only apply the former 
because we do not have a measure of the g — r color of 
the LMC. To apply their relationship, we need to convert 
our measured -F3.6, which is 2012 Jy for the the sum of 
all of the regions we studied, to L 3 . 6 in solar units. We 
use the zero point IRA C channel 1 calibration provided 
by IReach et all ff2005h (0 Ve ga magnitude corresponds 
to 280.9 J y) and a 3.6 /xm magnitude of the Sun of 3.24 
(jOh et all 120091 ) to calculate that L 3 . 6 for the LMC is 
3.7 x 10 9 L©. Applying their Eq. 2 results in an estimate 
of M* of 4.0 x 1O 1O M , which is a factor of eight larger 
than the dynamical mass. One potential source of the 
discrepancy is the choice of the Solar magnitude (because 
they do not quote their adopted value), although that 
seems unlikely to create an eight-fold difference. Another 
potential source is our measured L3.6, which could either 
be affected by errors or by contamination of luminous 
sources. We checked against errors by measuring the lu- 
minosity from a single aperture across the image (rather 
than summing the results from all of our subregions) and 
by comparing the lu minosity to the 3 . 6 ftm Tully-Fisher 
relation provided by iFreedman et ail ([201 ll ). Both com- 
parisons confirm that there is no significant problem with 
our measured luminos ity. The pr o blem may lie with ex- 
trapolating the fitted I Zhu et al~l (|2010D relation to the 
LMC, because adopting reasonable g — r colors for the 
LMC does decrease the discrepancy somewhat (although 
it does not bring down the estimated M* below Mb) 
or may reflect a metallicity dependence, which they ac- 
knowledge to have not addressed. 

Finally, we close with a discussion of the gap between 
Mo and M*, which is substantial (~ 3 x 10 9 Af Q or a 
factor of ~ 2.5). This "shortfall" could be ascribed to 
the ubiquitous dark matter, but it also leaves room for 
significant amounts of baryonic matter. We focus now 
on the impact of the choice of the IMF and implications 
for the study of unresolved stellar populations. There 
continues to be much debate over the appropriate func- 
tional form of the IMF and whether it is variable either as 
a function of galax y type (Ivan Dokkum fe Conrovll2"010l) 
or across redshift ([Pavel I2008D . There are two principal 
difficulties in applications of the stellar population mod- 
els in constraining the IMF in systems with unresolved 
stellar populations. First, the population of stars that 
dominate the total luminosity typically includes evolved 
stars, which are notoriously difficult to model. The mea- 
sured luminosity, relative to the model one, sets the nor- 
malization of the IMF. Second, the bulk of the inferred 



mass comes from low mass stars, which are not observed, 
leading to the stron g sensitivity of the implied ma ss on 
the form of the IMF. Ivan Dokkum fe Conrovl (|2010l ). us- 
ing a spectral signature of low mass stars, recently ar- 
gued (for giant elliptical galaxies) that the slope of the 
IMF might be significantly steeper than Salpeter, lead- 
ing to a vast reservoir of unappreciated stellar mass. Our 
study by-passes the first problem by counting stars rather 
than integrating luminosity, but is susceptible to the sec- 
ond one. For reference in c omparing the effe cts of IMF 
choice, an IMF favored by IBell et all (|2004f ) called the 
diet-Salpeter, because it reduces the total integrated stel- 
lar mass by limiting the number of low mass stars relative 
to Salpete r, results in 0.1 d ex less implied mass, while the 
use of the lChabrierl ([2003) mass function results in 0.25 
dex less implied mass. 

Our calculations above result in an inferred M^/L^^q = 
0.5 for the LMC (1.9 x lO 9 M /3.7 x 10 9 L Q ). This is al- 
ready a fairly low value of M^/L^ q, which would only 
decrease further if cither the diet-Salpeter or Chabrier 
IMFs were adopted. Given that the dynamical mass is 
already a factor of ~ 2.5 higher than the stellar mass, 
we contend that lowering the stellar mass even further 
is disfavored because we do not expect the dark mat- 
ter fraction to be signifi cantly larger than 60 % in the 
inner 3 kpc of the LMC (|Alves fc Nelson! 120001 ) . For ex- 
ample, adopting a Chabrier IMF (applying simple the 
0.25 dex correction rather than re-calcuating the SFH) 
results in the dark matter fraction increasing to ~ 80%. 
Such a high dark matter fraction is at the limit of the 
range found n ecessary to prod uce the baryonic Tully- 
Fisher relation (Mc Gau gh 2005) , and it is found in galax- 
ies of lower masses and surface brightnesses than the 
LMC. Going in the other direction, it is interesting to 
consider the botto m- heavy IMF suggested f or gia nt el- 
liptical galaxies by Ivan Dok kum fc Conrovl (|2010h . To 
fully consider the effects of an IMF that is different than 
Salpeter at all stellar masses, rather t han for stars with 
low m ass that are not constrained bv lHarris fc Zaritskvl 
(2009), requires recalculating the SFH with this new 
IMF. However, if we simply calculate the difference in 
stellar mass b etween Salpeter [x = — 2 .35) a nd the IMF 
suggested by Ivan Dokkum fc Conrovl IpOlCh (x = -3) 
when we normalize the two to have the same number 
of stars with M > 1M ), we find that the bottom- 
heavy IMF has a factor of 2.5 more mass (for stars with 
0.1 < M < 100 Mq). Given that our dynamical mass 
estimate is 5 x 1O 1O M0, we find that such a steep IMF 
is consistent with our measurements, but only if there is 
a negligible dark matter fraction over the inner 3 kpc of 
the LMC. 

In contrast, various other studies of late type galaxies 
contend that bottom heavy IMFs are inconsistent with 
dyn amical or lensing m ass constraints (see, for exam- 
ple, iBrewer et al.ll2012l). Given th e excellent data and 
analysis in the iBrewer et all (|2012D study, the only way 
we can envision reconciling the results is if the stellar 
population models are significantly underpredicting the 
luminosities used to match the observed luminosities. If 
the stellar models do underpredict ithe luminosity, that 
would in turn lead investigators to normalize the IMF's 
too high, resulting in a mass that violated dynamical 
constraints, unless the IMF is constructed to significantly 
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turn over at low masses. Intriguingly, some modelers are 
reaching analogous conclus ions for independent reasons 
(|Leitherer fc Ek strom 2 0111) . Alternatively, a systematic 
error lurks in our calculation of the stellar mass corre- 
sponding to the published SFH. Resolving this issue is 
manifestly paramount to many areas of astrophysics. 

4. SUMMARY 

We have used spatially resolved maps of stellar mass 
and IR flux in the Large Magellanic Cloud to 1) cali- 
brate a conversion between 3.6 and 4.5 /im fluxes and 
stellar mass, 2) examine potential approaches for using 
either 4.5 /im or 8 /an data to refine the estimates of 
stellar mass, and 3) provide empirical estimates of the 
uncertainty in such measurements. 

We find that: 

One can use measurements of the fluxes at 3.6 and 4.5 
/im, F3.6 and -F4.5, to estimate the corresponding stellar 
mass using 

Af — inS-65 T7i2.85 p-1.85 ( ® \ 

M.-10 F 3S F 4 . 5 ^—j 

where M* is in solar masses, F's are in Jy, and D is the 
distance to the source in Mpc. 

Although we were able to clearly identify deviations from 
the mean relationship between flux and stellar mass that 
correlate with stellar population variations, in particular 
with populations < 300 Myr old, and with hot dust, as 
traced by 8 fim emission, correcting for those sources of 
scatter, with data far superior to what will generally be 
available for most galaxies, results in marginal reductions 
in the scatter. Furthermore, we find that small fractions 
of such young populations (< 5%) are sufficient to sig- 
nificantly affect ^3.6- We conclude that it is difficult to 
reduce the scatter below what we find using broad, global 
measurements . 



The scatter, for our regions, which typically contain be- 
tween 10 5 - 5 and 10 6 5 M Q (and hence we presume compa- 
rable numbers of stars), is approximately 30% in mass. 
This scatter decreases as y/~N as regions with more stars 
are analyzed, resulting in an estimate of M* that has an 
internal precision of < a few percent for the entire LMC. 

Although the LMC contains a range of environments, 
covering strongly star-forming regions such as 30 Dor to 
quiescent regions, the relationship we provide may break 
down for even more strongly star-forming systems. We 
also have not explored the dependence of this calibration 
on metallicity. Within the LMC there is little variation 
in metallicity (Pagcl ct al. 1978), making it difficult to 
explore this issue within the context of this study. How- 
ever, one could extend this work to the SMC and to 
more distant galaxies for which resolved CMDs exist to 
establish the magnitude of the metallicity dependence. 
Finally, we have not directly addressed the uncertainty 
resulting from the adopted IMF. However, because this 
is a systematic uncertainty it will be more prominent 
in some uses of the relation, such as when total stellar 
masses are needed, than in other uses, such as when one 
is studying the distribution of stellar mass within one 
system. We have argued that our results can accommo- 
date the conjecture of a relatively bottom-heavy IMF for 
the LMC, such as Salpeter, and are in moderate con- 
flict with bottom-light IMFs such as the diet-Salpeter 
and Chabrier, which are usually advocated for late-type 
galaxies. Bearing these caveats in mind, we expect our 
calibrated conversion to be particularly useful in analyz- 
ing the large amount of extragalactic Spitzer and WISE 
data already in hand. 
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